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Abstract. A technique for constructing an infinite tower of pairs of "PT-symmetric 
Hamiltonians, H n and K„ (n — 2, 3, 4, . . .), that have exactly the same eigenvalues 
is described and illustrated by means of three examples (n = 2, 3, 4). The eigenvalue 
problem for the first Hamiltonian H n of the pair must be posed in the complex domain, 
so its eigenfunctions satisfy a complex differential equation and fulfill homogeneous 
boundary conditions in Stokes' wedges in the complex plane. The eigenfunctions of 
the second Hamiltonian K n of the pair obey a real differential equation and satisfy 
boundary conditions on the real axis. This equivalence constitutes a proof that 
the eigenvalues of both Hamiltonians are real. Although the eigenvalue differential 
equation associated with K n is real, the Hamiltonian K n exhibits quantum anomalies 
(terms proportional to powers of h). These anomalies are remnants of the complex 
nature of the equivalent Hamiltonian H n . For the cases n = 2, 3, 4 in the classical limit 
in which the anomaly terms in K n are discarded, the pair of Hamiltonians H Ui classical 
and K nt classical have closed classical orbits whose periods are identical. 
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1. Introduction 



This paper discusses the infinite class of higher-derivative Hamiltonians 

ff„ = f"-7(«) ( " 2) (n = 2, 3, 4, ...), (1) 

where rj and 7 are arbitrary positive real parameters. These Hamiltonians are all VT 
symmetric; that is, they are symmetric under combined space reflection V and time 
reversal T . The literature on VT symmetry is extensive and growing rapidly; some 
early references [U [21 [3l H] and recent review articles EJ [7] provide background for 
this work. 
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The Hamiltonians in ([T]) are special because they are spectrally identical to another 
class of Hamiltonians K n , which have entirely real spectra. For example, for the case 
n = 2 the Hamiltonian 

B^^f-^ (2) 

has a real positive discrete spectrum and this spectrum is identical to the spectrum of 
the conventionally Dirac Hermitian Hamiltonian 

k 2 = ^ + A^ + h^p. (3) 

The spectral equivalence of these two Hamiltonians has long been known and has been 
examined in many papers El [101 El CGI CCU E] . 

Spectral equivalences between pairs of non-Hermitian and Hermitian Hamiltonians 
have been discussed in the past [151 EH] an d a number of approximately equivalent pairs 
of Hamiltonians have been constructed (see, for example, Refs. [T7l [TH]). A few other 
exactly equivalent pairs of Hamiltonians have been found [191 EH EH E2, E31 EH EH] . 

A brief demonstration that the Hamiltonians in (T5]) and ([3]) are spectrally identical 
is given in Ref. |llj . The prescription introduced in Ref. [H] makes use of several 
elementary transformations of differential equations. We review this analysis here. 

To begin, we construct the coordinate-space eigenvalue problem for H2 in ([2]) by 
substituting p = —ih-^. The formal eigenvalue problem Hztp = Eip then becomes the 
differential-equation eigenvalue problem 

h 2 

i)"(x) - jx 4 i[)(x) = EMx). (4) 

2m 

The eigenfunction ip(x) is required to satisfy homogeneous boundary conditions inside 
a pair of Stokes' wedges in the complex- 2 plane. These wedges lie below and adjacent 
to the real-x axis and have angular opening tt/3, as shown in Fig. [TJ 

To establish that the eigenvalues of (jlj) are also eigenvalues of Ki in 0, we follow 
the simple four-step recipe introduced in Ref. [TT] : 

Step 1: Map the integration path onto the real axis by the change of variable 



x = —1i\J\ + it, t E (— 00,00). (5) 

As t runs along the real axis from —00 to +00, x runs along the complex contour that 
asymptotes in the Stokes' wedges shown in Fig. [U Under the change of independent 
variable in ([5]), the derivatives transform as follows: 



a =VT+Tt± (6) 

Thus, with the substitution if)(x) = (pit), (jlj) becomes 

- £ [a + itw{t) + \i<j>'{t)\ - 167(1 + it?m = EM). (8) 

Step 2: Take the Fourier transform of jSJ). The Fourier transform of <f>(t) is 

f(p)= [ dte^it), (9) 
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so the Fourier transform of (|S]) is achieved by making the replacements 

4 and 



if-ft± and 4--$, (10) 



and we obtain 

- I6 7 n 2 /"(P) + A (£ - 32 7 ) /'(p) + + g _ ie 7 ) /(p) 

= S/(p). (11) 
Step 3: Introduce the new dependent variable g(p) by 

f(p) = Q(p)g(p) (12) 

and choose Q(p) to remove the one-derivative term and thereby change (lllj) to 
Schrodinger form. The condition on Q(p) that eliminates terms containing g'(p) is 

Q'(p) _ p 2 _ l f 1 q\ 

Q(p) 64m7ft ft" 

Differentiating fflBT) . we get 

Q(p) 32m-yh ~ \64m"/h hj v ^ 

and substituting ^2l - [H|) into (ITTj) . we obtain the Schrodinger equation 

- 167»V(P) + + il) *(p) = (15) 
Step 4: Finally, rescale the independent variable p: 

p — ► A/32m 7 p. (16) 



Im x 




Figure 1. Stokes' wedges in which the eigenfunctions f° r the — x 4 eigenvalue 

problem (TJJ are required to vanish as |x| — > oo. These wedges border on, but do 
not include the real axis. A possible path along which one can solve the eigenvalue 
differential equation is shown. This path must asymptote inside the Stokes' wedges. 



Exact Isospectral Pairs of VT -Symmetric Hamiltonians 



4 



Im x 




-1.5 -1.0 -0.5 



0.5 1.0 1.5 



Rex 




Figure 2. Classical "PT-symmetric trajectories in the complex- x plane for the 
Hamiltonian -H2, classical with m = 1/2 and 7=1. For these trajectories the energy 
E = 1 and the trajectories start at 2i and — 2i. The period of the motion is 
T = 1.854 074 68... . 



This rescaling makes (j!5p resemble the original Schrodinger equation in (J3J): 

- £</"(p) + (47P 4 + 9(P) = Eg{p). (17) 

The Schrodinger equation f|T7|) is associated with the Hamiltonian in ([3]). Since the 
differential equation ( TITl) is real and the boundary conditions on the eigenfunctions g(p) 
are imposed on the real axis, it follows that the eigenvalues are real. This constitutes a 
rigorous proof that the non-Hermitian Hamiltonian H2 in ([2]) has a real spectrum. 

The term in ([3]) that is proportional to h is a quantum anomaly that arises as a 
remnant of the complex boundary conditions on the eigenfunctions associated with (j4]). 
(These boundary conditions are illustrated in Fig. [TJ) In the classical limit h — > 0, we 
obtain the classical Hamiltonian 

K 2) classical = ^ + 4 7 p 4 . (18) 

It is an interesting but previously unnoticed fact that at the classical level the two 
Hamiltonians ([2]) and ffTBl are equivalent. To demonstrate this equivalence we show that 
the classical trajectories determined by these Hamiltonians in complex-coordinate space 
have identical periods. We plot in Fig. [2] two classical orbits of if 2 , classical for m = 1/2, 
7 = 1, and energy E — 1. Note that there are four turning points, which are located at 
x = e ln ^, x = e 3 * 71 "/ 4 , x = e 5 * 71 / 4 , and x = e 7 * 71 "/ 4 . A trajectory starting at x = 2i gives a 
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Figure 3. Classical "PT-symmetric trajectories in the complex-x plane for the 
K% t classical Hamiltonian in p8[) . The period of this orbit is identical to the periods 
of the orbits shown in Fig. 



closed orbit in the upper-half plane, and a trajectory starting at x — —2i gives a closed 
orbit in the lower-half plane. Both orbits have the same period T, where 

= v / 2^r(i/4) = 854 074 68 (19) 

4r(3/4) v ; 

We calculate T by solving Hamilton's equations x = = 2p and p = — ^ = Ax 3 . 
Eliminating p, we express the period T as the contour integral 



2V1 + 



x 



We then use Cauchy's theorem to distort the contour to one that goes along a ray 
from the origin to a turning point, encircles the turning point, and follows the ray back 
to the origin. The contour then continues along another ray out to the other turning 
point, encircles this turning point, and follows that ray back to the origin. The resulting 
integral becomes a standard representation for a Beta function. 

In Fig. [3] we plot the classical orbit for a particle of energy E = 1 described by 
the Hamiltonian K 2> classical in ([18]) . Note that there are now two and not four turning 
points located at x — ±1. The figure shows a closed periodic classical orbit that begins 
at x = 2i. The period T of this orbit is exactly that given in (fl9|) . 

Figure [3] exhibits a remarkable new topological feature that is not found in the 
complex classical trajectories of Hamiltonians of the form p 2 + V(x), namely, that the 
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Figure 4. Behavior of a classical trajectory as it approaches a turning point of a 
Hamiltonian of the form p n + V(x) for n = 2, n = 3, and n = 4. When n = 2 the 
trajectory executes a 180° U-turn, when n = 3 the trajectory makes a 240° turn, and 
when n — 4 the trajectory makes a 270° turn. The wiggly lines indicate that the 
turning point is a branch point of an n-sheeted Riemann surface. Turning points of 
the n = 3 type are shown in Figs. [S] and [T^] and turning points of the n = 4 type are 
shown in Figs. [3] and [5] 

classical orbit makes a 270° loop around the turning points. There have been many 
studies of complex classical systems [2] [26], [271 [23 [29], [30] [31] and in previous numerical 
studies of complex trajectories the Hamiltonians that were examined were quadratic 
in the momentum p. When the momentum term in the Hamiltonian is quadratic, the 
trajectory always makes a 180° U-turn about the turning points. In contrast, with 
Hamiltonians of the form p 3 + V(x), the classical trajectory makes a 240° turn and for 
Hamiltonians of the form p A + V(x), the classical trajectory makes a 270° turn about 
the turning points. These behaviors are illustrated in Fig. [H 

To understand the angular rotation of a classical trajectory about a turning point, 
let us consider a Hamiltonian of the form H = p n + V(x). From the first of Hamilton's 
equations, x — ^ = np n_1 , we can eliminate p from the Hamiltonian and we find that 
a particle of energy E obeys the equation x = n[E — V(x)]( n ~ 1 ^ n . We can now see 
that the turning point is also a branch point of an n-sheeted Riemann surface. Note 
that as n — >• oo, a classical trajectory that approaches a turning point will make a full 
360° turn and continue going in the same direction! (We encounter many examples of 
such higher-order turning points throughout this paper. However, we emphasize that, 
in general, the properties of higher-order classical turning points are not easy to predict 
analytically, and require the kind of sophisticated analysis used in catastrophe theory.) 

The main objective of this paper is to show how to construct isospectral pairs of 
Hamiltonians, one of which is associated with a complex eigenvalue problem while the 
other is associated with a real eigenvalue problem. There have been many attempts 
to find isospectral pairs of Hamiltonians of the form p 2 + V(x) by using the four-step 
differential-equation procedure outlined above. However, all such attempts have proved 
fruitless. In this paper we consider a wider class of Hamiltonians in which p 2 is replaced 
by p n , and we suggest by means of illustrative examples (rather than by presenting a 
proof) that for each integer value of n it is now possible to construct such an isospectral 
pair of PT-symmetric Hamiltonians. (To construct a proof, one would have to follow 
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the procedure detailed in the illustrative examples in Sees. [3] and 01) The eigenfunctions 
for the first member of the pair H n satisfy boundary conditions in Stokes' wedges in 
the complex plane. However, the eigenfunctions for the second member of the pair 
K n satisfy a real differential equation with homogeneous boundary conditions given 
on the real axis. Therefore, our construction constitutes a rigorous proof that entire 
eigenspectrum of the complex Hamiltonian H n is real. 

The Hamiltonian K n that is spectrally equivalent to H n has quantum anomaly 
terms containing powers of h up to n — 1. If we discard these anomaly terms, we obtain 
a pair of classical Hamiltonians that are equivalent in the following sense: For each 
closed periodic classical trajectory of the first Hamiltonian, it appears that there exists 
a closed trajectory of the second Hamiltonian that has exactly the same period. 

This paper is organized as follows: In Sec. [2] we present some general introductory 
calculations. Then, in Sees. [3] and H] we treat the cases n = 3 and n = 4. Finally, in 
Sec. [5] we make some concluding remarks. 



2. General treatment 

In order to carry out Step 1 on H n in ([3]) we must generalize the mapping in (Q: 

x = --(l + it) a , te (-00,00), (21) 
a 

where a = 1/n. Note that this change of variable reduces to ([5]) when n = 2 (a = |). 
Now, as t runs from —00 to +00 along the real axis, x runs along a complex contour 
that is appropriate for the Hamiltonian H n . 

The transformations from x derivatives to t derivatives can then be written as 

± ={1 + it y-«A (22 ) 
(A) 2 = (1 + it)*-** (A) 2 + i(1 _ Q)(1 + it f-^ d ( (23) 

(£) 3 = (1 + it)*-** (A) 3 + 3,(1 -«)(! + ,t) 2 "- (A) 2 

-(l- a )(l-2a)(l + tt) 1 - 3a ± (24) 

(£) 4 =(l + ^) 4 - 4Q (^) 4 + 6,(l-a)(l + ,t) 3 --(A) 3 
-(l-«)(7-ll«)(l + ^) 2 - 4a (A) 2 

- i(l - a)(l - 2a)(l - 3a)(l + it) l ~ ia f t . (25) 

This set of equations generalizes (JBj) and ([7j); fl22l) and (1231) reduce to (jSJ) and (j7]) when 
a = |. In the next two sections we show how use these results to transform a general 
class of Hamiltonians of the form ([1]). Note that the Hamiltonian in ([2]) corresponds to 
n = 2. We begin with the simplest generalization, namely, n — 3. 

3. Case n = 3 

In this section we consider the Hamiltonian (CQ) for the case n = 3: 

H 3 = np 3 - i-fx 9 . (26) 
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Figure 5. Stokes' wedges for the solutions to the differential equation (|27|) . The 
boundary conditions on the solution to this eigenvalue problem require that ip(x) — > 
as \x\ — ► oo with arg£ inside the two shaded wedges in the lower-half plane. Two 
integration contours inside these wedges are shown, one of which runs from an angle 
of — 57r/6 to an angle of — ir/Q and the other from — 2n/3 to — n/3. When we set a = | 
in (I2ip and allow t to run from — oo to oo, the variable x follows the latter curve in 
the complex- x plane. 



The corresponding differential-equation eigenvalue problem H^ijj = Ei/j is 

iH 3 T] {£) 3 ^(x) - i-yx 9 ip(x) = Eip(x). (27) 

To solve this boundary-value problem the integration path must lie in Stokes' 
wedges. For the above equation these wedges do not include the real axis. The 
wedges are determined by the asymptotic solutions to (!27l) . which for large |x| have 
the form exp(cu;a; 4 ), where c = ^(7/??) 1//3 is a positive constant and uj is a cube root 
of unity: u 3 = 1. From this asymptotic behavior we see that each of the Stokes' 
wedges has angular opening it/ 4. We impose boundary conditions in the wedges 
— 7r/8 > argx > — 3n/8 and — 5n/8 > argx > —7n/8. (These wedges are the shaded 
regions in the lower-half plane shown in Fig. 0) Requiring that ip(x) vanish as \x\ — > oo 
with arg x in the above two wedges eliminates two of the three possible solutions to (1271) 
and keeps only the solution whose asymptotic behavior is given by exp(er 4 ). 
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We now follow Step 1 of the procedure explained in Sec. [T] and make the change of 
variable in (121"]) with a. — |. From (124"]) we then obtain the differential equation 

ih 3 v [(i+.t) 2 (^) 3 +2i(i+it) ay 

- 3 9 7 (l + it) 3 0(t) =E<j>(t). (28) 

This differential equation is to be solved on the real axis in the t variable and the solution 
satisfies homogeneous boundary conditions on the real axis. Note that as t runs from 
— oo to oo, the variable x runs from an angle of — 27r/3 to an angle of — tt/3 in the 
complex- x plane. This integration in the complex- x plane is shown in Fig. 

Following Step 2, we perform the Fourier transform in ([9]) and use ( TTUj) to obtain 

2 



/(P) 



- 3 9 7 (l + ^) 3 /(p)=^/(p)- (29) 

We then simplify f j29|) by transforming it to a new differential equation that does 
not have a second-derivative term. To do so we let f(p) = Q(p)g(p), as in Step 3 of 
Sec. [TJ The condition on Q(p) that eliminates the second-derivative term is 

~M - _I _ jzeL. (30) 



Q(p) ft 3 10 7ft 

and upon differentiation we get 



Q"(p) 
Q(p) 

Q"'(p) 
Qip) 



j_ 

r, 2 



VP 



+ 



3 9 7?i ^ 3 20 7 2 ^ 2 ^ 3 10 7^ 2 ' 



+ 



2rjp s 



W 3^E ^ 3 18 7 2 ^ 2 ^ 



ft2 3 19 7 2 ^3 



VP 

WW 



3 30 7 3 /i 3 • 



(31) 



The resulting equation for g(p) is 
- 3 9 7 n^"'(p) + 

+ 



4q 2 hp 5 
3 io 7 



2rih 2 p 
9 



77/iV) 0'(p) 



2??- i p 9 
3 21 7 2 



This completes Step 3. 

Finally, we perform the scaling 



P 



277 1 / 3 



t^rP 



(32) 



(33) 



and obtain 



- V h 3 g"\p) + (-27r7 2 /y/ 3 fry + 243r / 1 / 3 7 2 / 3 ^ 6 ) </(p) 

+ (972r/ 1/3 7 2/3 ^ 5 - 6r/ 2/3 7 1/3 n 2 p + I4587/) = £#(p). (34) 

This completes Step 4 and from this result we make the ansatz j- — > |s to identify the 
Hamiltonian that is equivalent to that in (1261) : 

K 3 = ir]x 3 + i (-277? 2/3 7 1/3 frp 2 + 243r] 1/3 7 2/3 p 6 ) x 

+ (972r/ 1/3 7 2/3 ^ 5 - 6r/ 2/3 7 1/3 n 2 p + 14587p 9 ) . (35) 

Observe that there is an h 2 anomaly term as well as two H anomaly terms. Therefore, 
in the classical limit h — > 0, K 3 in ( 135]) becomes 

K 3 , classical = %^ + 243w7 1/3 7 2/3 A + 1458 7 p 9 . (36) 



Exact Isospectral Pairs of VT -Symmetric Hamiltonians 



10 



Im x 
1.5^ 




Figure 6. A classical "PT-symmetric periodic trajectory in the complex- a; plane for 
the Hamiltonian classical = P 3 — i% 9 for a particle of energy E — 1. The trajectory 
begins at x = 0.5 i and proceeds to wind around six of the nine turning points. It never 
crosses itself because each of the turning points is a cube-root branch point, and thus 
the trajectory visits the three sheets of the Riemann surface, which are numbered 0, 
1, and 2. The period of the motion is T = 3.377 981 999 015 . . . . 

We now argue that the two classical Hamiltonians, H 3> classical and K 3> classical) are 
equivalent by computing numerically the classical trajectories for each Hamiltonian and 
verifying that trajectories having the same energy have the same periods. In Fig. [6] we 
plot a classical trajectory for H 3 c i aS sicai with 77 = 1 and 7 = 1 for a particle of energy 
E = 1 starting at x — 0.5 i. 

The trajectory in Fig. [6] has an elaborate structure. It winds around six of the nine 
turning points in the negative (clockwise) direction. Because each of the nine turning 
points is also a branch point of cubic type, with each wind the direction of the path 
changes by 240°, as shown in Fig. |H Thus, after the path leaves its starting point at 
x = 0.5 i, it winds in the negative direction and crosses from sheet zero to sheet two 
of the Riemann surface (which is equivalent to sheet —1). The number of the sheet is 
indicated in a circle on Fig. [61 The path then heads upward and to the left, but it does 
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not cross itself because it is on the second sheet and no longer on the zeroth sheet. The 
path then proceeds to visit each sheet of the Riemann surface twice before returning to 
its starting point at x = 0.5 i on the zeroth sheet. 

We can calculate the period of this orbit exactly by distorting the path shown in 
Fig. [6] into a path that begins at the origin, travels outward to a turning point along a 
ray, encircles the turning point, and then returns to the origin, and then repeats this trip 
five more times. However, in evaluating the complex line integrals we must be careful to 
remember that the phase of the integrand changes as the path encircles the turning point. 
Upon adding together the contributions from all six turning points we obtain the exact 
result T = 2[cos(vr/18) + cos(77r/18)]r(10/9)r(l/3)/r(4/9) = 3.377981 999 015 .. . . 

In Fig. [TJwe plot the classical trajectory of a particle of energy E = 1 governed by 
the Hamiltonian K% t classical with r\ = 1 and 7 = 1. This trajectory begins at x = 0. The 
time for the particle to follow the path shown in this figure is T = 1.125 994 . . ., which 
is precisely 1/3 of the period of the closed orbit shown in Fig. [61 We emphasize that 
the trajectory of this orbit is not yet closed. Figure [7] shows that the particle must trace 
this path three times before the orbit closes, and the time required to do this is exactly 
the period of the orbit shown in Fig. [6j 

One might think (wrongly!) that there should be a total of six rather than three 
turning points on Fig. We determine the positions of the turning points from the first 
of Hamilton's equations 

x = ^-K 3i classical = 3 6 2tp 5 x + 3 8 2p 8 . (37) 
op 

The turning points are the points where x — 0, and this seems to occur when p = 
and when ix + 9p 3 = 0. However, if we then substitute each of these conditions into 
the equation K 3t classical = 1, we find that only the first of these conditions is consistent, 
so we learn that there are three turning points situated at the three roots of —i. The 
second condition is inconsistent. 

4. Case n = 4 

In this section we consider the Hamiltonian 

H 4 = - 7 £ 16 . (38) 

The substitution p = —ih-^ converts the formal eigenvalue problem H^ip = Eip into the 
differential-equation eigenvalue problem 

(£) 4 ^(x) - 7x lf V(x) = Eip(x). (39) 

The asymptotic behavior of the eigenfunctions ip{x) has the form exp(icr 5 /5), where 
c = (7/77) 1//4 is a positive constant. Thus, the eigenfunctions vanish exponentially fast 
in a pair of PT-symmetric wedges in the lower-half complex-x plane centered about 
-37T/10 and -7vr/10. 
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Figure 7. A classical 'PT-synimetric trajectory in the complex- x plane for the 
Hamiltonian K^ t classical in (|35[) with rj = 1 and 7=1. The trajectory shown is that of a 
particle of energy E = 1 beginning at x — —i. The trajectory follows an elaborate and 
complicated path as it winds around the three turning points, which are indicated by 
dots. The inset shows a blown-up version of the trajectory near the origin. Note that 
the trajectory never crosses itself because each of the turning points is also a cube-root 
branch point, and thus the trajectory lies in a three-sheeted Riemann surface. The 
period of the orbit is T — 1.125 994 . . ., which is exactly 1/3 of the period of the closed 
orbit shown in Fig. [6] It is crucial to observe that the trajectory shown is not closed] 
the inset indicates that if the trajectory is on sheet just to the right of the origin, it 
returns to this point on sheet 1. Thus, it must trace this path two more times before 
it closes. When it finally closes, the period of its orbit agrees exactly with that of the 
trajectory shown in Fig. [6] 



A, A 

32 dt 



<p(t) 



Step 1: We map this differential-equation eigenvalue problem onto the real axis 
by taking a = | in (pTj) so that x = (1 — it) 1 / 4 . The mapped eigenvalue problem is 

rjh 4 [(1 + rtf (A) 4 + + itf (A) 3 _ §(1 + it) ( d y _ 

-A 16 1 {l + it) 4 (t ) {t)=E(P{t). 

Step 2: Taking the Fourier Transform of fj40|) . we get 

3 „ / . \ 2 



(40) 



r]h 4 



^ + ' l dp J W 2\ 1 ^ a dp J fl 3 16 y ' l dp J h 2 32 h 

-4 16 7 (l + h±yf(p) = Ef(p). 



f(p) 



(41) 



Step 3: We transform the dependent variable using f(p) = Q(p)g(p) and choose 
Q(p) to eliminate the third-derivative term g"'(p). The function Q(p) that does this is 



Q(P) = exp (-1 + J^) ■ 



(42) 
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The result in (142 j) simplifies (|4ip to 

- 4 16 /z 4 7 </» + 3 -^g(p) + [8 P 12 g'(p) + Mp^g 



+ [2Ap s g"(p) + 2A0p 7 g'(p) + A83p s g(p)} 

- § i-48p 3 g"(p) - Qp 2 g'(p) + 87pg(p)} = Eg(p). (43) 
Step 4: Next, we perform the scaling transformation 

P -> ^?7^P (44) 
and obtain the real eigenvalue problem 

- vtfg""(p) + 4 12 3 7 p 16 £Go) + M 9 7 3/ V /4 [8p 12 g'(p) + 54p n (7(p)] 
+ h 2 4 5 ^q [24p 8 g"(p) + 240p 7 g'(p) + 483p 6 g(p)] 

- 8^ 3 7 1/4 7] 3/4 [-48p 3 ^"(p) - 6pV(p) + 87pg{p)] = Eg(p), (45) 

whose boundary conditions are posed on the real axis. Finally, we identify the 
Hamiltonian that gives rise to this eigenvalue problem by substituting ^ — > jx: 

K A = - r,x A + 4 12 3 7 p 16 + 4 9 7 3/ V /4 (8^ 12 x + hAhp 11 ) 
+ 4 5 v ^ (-24j3 8 x 2 + 240hip 7 x + 483h 2 p 6 ) 

- 87 v y /4 (A8hp 3 x 2 - 6ih 2 p 2 x + 87h 3 p) . (46) 

We have thus shown that the eigenvalues of K A are real and are identical to the 
eigenvalues of H&. 

Notice that K± has first-, second-, and third-order anomaly terms. In the classical 
limit S^Owe obtain 

^.classical = -2 13 3V^A 2 - V x* + 2 2 H 1 ^ly 2 x + 2 24 3 7 p 16 . (47) 

We will now demonstrate that the two classical Hamiltonians, classical and classical) 
are equivalent by showing that the classical periods of the motion are identical. 

Let us first examine the complex particle motion due to classical for the case 
7] = 1 and 7=1. In Fig. [8] we plot a closed orbit of a particle of energy E = 1 
that begins at 0.8i. The orbit takes the approximate shape of a rectangle and makes 
270° turns at the four corners. The period of this orbit is T = 1.697 723 019 533. . . . 
We can calculate this period exactly by distorting it into rays that go from the 
origin out to each of the turning points and back to the origin, taking great care 
to keep track of the phase of the integrand. The exact formula for the period is 
T = [cos(3vr/16) + cos(5vr/16)]r(l/4)r(17/16)/r(5/16). 

Next, we examine the classical trajectories of -^classical for the case r\ = 1 and 
7 = 1. In Fig. [9] we plot the classical trajectory for a particle of E = 1 that begins at 
x — — i. There are four turning points of fourth- root type. However, in this figure it is 
not possible to see the detailed structure of this trajectory near the origin. Therefore, 
in Fig. [10] we plot a blow-up of the square region in Fig. [9] near the origin. 

In Fig. [10] we can now see the four turning points and a rather complicated 
trajectory. The regions surrounding the upper pair of turning points have a highly 
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Figure 8. A classical "PT-symmetric periodic trajectory in the complex- a; plane for 



the Hamiltonian classical = P 



,.16 



for a particle of energy E = 1. The trajectory 



begins at x — 0.8i and proceeds to wind around four of the 16 turning points. It never 
crosses itself because each of the turning points is a fourth-root branch point, and thus 
the trajectory visits the four sheets of the Riemann surface, which are numbered 0, 1, 
2, and 3. The period of the motion is T = 1.697 723 019 533 ... . 



detailed structure, and thus the square region surrounding the upper-right turning point 
is blown up again in Fig. [11] (left side). The vicinity of the turning point contains 
additional structure, and this region must be blown up still more in Fig. CD] (right side). 

5. Summary 

We have shown by using a number of examples that it is possible to construct an infinite 
tower of pairs of isospectral Hamiltonians H n and K n , for which the first member of 
the pair is a complex "PT-symmetric Hamiltonian. The differential-equation eigenvalue 
problem for the second Hamiltonian is entirely real, and therefore the eigenvalues of 
both Hamiltonians are real. The second member of the pair has quantum anomalies of 
order ft n_1 . We have also shown that at the classical level, where the anomaly terms 
in K n are discarded, the Hamiltonians are equivalent by demonstrating that they have 
closed orbits of the same period. 
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Im x 




Figure 9. A classical trajectory in the complex- x plane for the Hamiltonian Ki t classical 
in (147|) with 77 = 1 and 7 = 1. The trajectory represents a particle of energy E = 1 that 
begins at x — e~™/ 4 . The trajectory winds around the turning points in a complicated 
fashion (see nested insets in Figs. ITUl and ITTj) but it never crosses itself because it is on 
a multisheeted Riemann surface with each of the turning points being a fourth-root 
branch point. The trajectory returns to its starting point in time T = 0.848 86, but it 
is not a closed trajectory; the trajectory begins on sheet and returns on sheet two of 
the Riemann surface. To complete its periodic motion the trajectory must follow the 
same path one more time. The time needed to complete this double loop is exactly 
equal to the period of the orbit shown in Fig. [8l 

Many unsolved problems remain. For example, we do not know if for every 
classical orbit of the Hamiltonians -ff n , classical there is a corresponding classical orbit 
of the Hamiltonians K n> classical of exactly the same period. Because complex coordinate 
space is a Riemann surface having an elaborate sheet structure, it is possible to find 
various classical orbits having many different periods. For example, for the classical 
Hamiltonian -£^3, classical with rj = 1 and 7 = 1, we find that a particle of energy 
E — 1 that starts at x — 1 has a closed trajectory of period T = 6[cos(7r/18) + 
cos(57r/18)]r(10/9)r(4/3)/r(4/9) = 4.143 708 353 ... . This orbit is shown in Fig. [12 

We calculate the period of this orbit exactly by distorting the path shown in Fig. [12] 
into a path that begins at the origin, travels outward to a turning point along a ray, 
encircles the turning point, and then returns to the origin, and then repeats this trip 
five more times. In the complex line integrals the phase for each path advances as the 
path encircles the turning point. Upon adding together the contributions from all six 
turning points we obtain an analytic expression for the period of the motion. Although 
we are convinced that one exists, we have not yet been able to find a corresponding 
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Figure 10. A blown-up view of the trajectory shown in Fig. [H] (see inset in Fig. [5]). 
The four fourth-root turning points are visible in this figure. However, the trajectory 
has a complicated structure near the upper two turning points, and Fig. [Tl] shows a 
detail of the small region surrounding the upper-right turning point. 




Figure 11. Left: A blown-up view of a portion of Fig. fTOl (see inset in Fig. [T0|) . Right: 
A blown-up view of the inset on the left side. 
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Figure 12. A classical "PT-symmetric periodic trajectory in the complex- a; plane for 
the Hamiltonian -ff 3, classical — p 3 — ix 9 for a particle of energy E — 1. The trajectory 
begins at x — 1 and proceeds to wind around six of the nine turning points. It never 
crosses itself because each of the turning points is a cube-root branch point, and thus 
the trajectory visits the three sheets of a Riemann surface, which are numbered 0, 1, 
and 2. The period of the motion is T = 4.143 708 353 ... . 

orbit of the Hamiltonian classical- 
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